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Nowadays, random telegraph signals play an important role in integrated circuit 
performance variability, leading for instance to failures in memory circuits. This 
problem is related to the successive captures and emissions of electrons at the many 
traps stochastically distributed at the silicon-oxide (Si-S;02) interface of MOS tran- 
sistors. In this paper we propose a novel analytical and numerical approach to 
statistically describe the fluctuations of current due to random telegraph signal in 
time domain. Our results include two distinct situations: when the density of inter- 
face trap density is uniform in energy, and when it is an u-shape curve as prescribed 
in literature, here described as simple quadratic function. We establish formulas for 
relative error as function of the parameters related to capture and emission proba- 
bilities. For a complete analysis experimental u-shape curves are used and compared 
with the theoretical aproach. 

Low frequency (LF) noise is a performance limiting factor for deep sub-micron CMOS 
devices ^]. In these devices, LF noise is dominated by multiple Random Telegraph Signals 
(RTS) (see jj] for a discussion on RTS in S;-Si0 2 interfaces). The MOS transistor employed in 
analog and digital designs drives a current Ids which is a function of the potential differences 
Vr>s and Vqs applied respectively between drain and source and gate and sourcejlj]. RTS 
noise is due to succeeding electron capture and emission by a number of N tr traps distributed 
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according to a Poisson distribution at the Sj-Sj0 2 interface, as it can be observed in plot [T] 
(a). This phenomena causes oscillations in the transistor current Ipg over time, even with 
Vds and Vqs constant. A typical variation in the current caused by one single trap is shown 
in figure [T] (b). 




FIG. 1: Plot a: Physical mechanism for the generation of LF-noise; Plot b: This figure shows the 
effects caused by only one trap according to physical mechanism described in plot a. A discontinuity 
in the electrical current corresponds to the absorption or emission of an electron. 

Noise performance may strongly vary between different devices on a same chip, and even 



between different operation points of a single transistor. Recent results {9] claim these effects 
responsible for fails in flash memories. For future technologies, RTS may become relevant 
also for the performance of combinational digital gates, as inverters and other logic gates. 
The propagation delay of a gate depends on the capability of its transistors to drive current. 
A smaller current driven by the transistor means larger propagation delay, which may lead to 
timing violations (failures). Hence, the variation in transistor drive current due to RTS may 
lead to circuit failures in future technology generations, and statistical modeling of random 
telegraph signal is demanded. But until now statistical models for RTS focused on frequency 
domain, for instance refer to j^] for a statistical RTS model in frequency domain. This is 
suitable for analog circuits, whose design and analysis is performed in frequency domain. 
However, for digital circuits an appropriate time domain statistical analysis is needed, since 
these circuits are analyzed and designed using time domain metrics. So, as RTS is becoming 
a major issue on upcoming digital integrated circuits, an approppriate RTS noise model for 
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time domain is mandatory for the correct modeling of the behavior of fugure digital circuits. 

Aiming to address this issue, this work presents for the first time in literature a compre- 
hensive model for RTS in time domain. We derive analytical expressions for the average and 
relative error of the transistor current. We also analyze the joint contribution of all the traps 
to the transistor current, instead of analyzing only the current fluctuation due to one single 
trap, which is more suitable for analyzing the total current fluctuation of the transistor. Our 
model considers the density of states in the band-gap analyzing three possible situations: 

1. The density of states is uniform in the band-gap range. This is based on previous 
results from literature (see for example 

2. The density of states curve measured from experimental data, which was published 



in jg] and shows that the density of states follows an u-shape curve. We performed a 
polynomial fit with a 8-degree polynomial which fits very well with the experimental 
data. 

3. The density of states modeled as a quadratic function in the band-gap range. This is 



an approximation for the experimental data of js], which is an u-shape curve. 

In the random telegraph signals, cr, = 0, 1 denotes the state of the i-th trap (0 =empty 
or 1 =occupied), the Fermi-Dirac statistics governs the probability of transition: 

_ , n * . dt dt 

Pr ^ = ->■ <7i = l)dt = = 

10P- [1 + exp(-gi)J T W 

dt dt 

PrlCTj = 1 — > Oi — 0)dt = ; : — — = —pr 

1 ' 10w [1 + expfe)] r « 

where tOP and Te^ are respectively the average time of capture and emission and = 
(Ej) —jj) I /cbT where Ej) is the energy ( within the band-gap) of the i— th trap, /x is the fermi- 
level energy, ks = 1.3806568 x 10~ 23 JA;^ 1 the Boltzmann constant and T is temperature. 
In a first model — Q < qi < Q is considered a uniform random variable and Tc \ are 
identically distributed, i.e., (rj ; ) = (r c ) and (re^) = (r e ) for i = 1,2, ...,N t 
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Here, Pi is also a random uniform variable within an interval p m [ Q < Pi < p max and in 
this case in the frequency domain, we can establish an important connection. It is not 
difficult to show (see for example [4H3]) that the power spectrum density correponding to 
the noise from the i-th trap is a Lorentzian function Si(fi) = (Aj/fi) [1 + (/ / fi) 2 } 1 where 
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fi = + is the corner frequency corresponding to the trap and Ai is its amplitude. 

From that we conclude fi = 10~ Pl and due to this fi is uniformly distributed in a log 10 
scale resulting in a probability distribution h(fi) = [In 10 (p max — Pmin) fi] 1 for the corner 
frequencies (this assumption will be used from now on in this paper). 
From this approach, we can calculate 

Prfo = l) = + 1) 



r c + r e 

and thus the average current due to a single trap is computed by 



h = Si ■ Prfo = 1) + • Prfa = 0) 



The amplitudes <5j are also random variables and our results will be dependent on its first 
and second moments, respectively (5) and (S 2 ). The next step is to consider the contribution 
of all traps, because by doing so one can compute the fluctuation of the current passing 
through the channel of the MOS transistor. The average current, taking in consideration all 
random sources, is written as: 

N tr 



1 / Ntr,T c ,T e 



S ( £ N tr 



N N tr( 



N tr =0 



tr- 



I 



1+e-i 
2(l+cosh q) 



dq 



N (S) 
2Q 



W [2Q + lne-«] 
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where N = c (j5 max — p min ) and c = In 10 Nd ec W L, where W and L are the device dimensions 
and p = In 10 N dec is the density of traps by area unit and by decade frequency in log-scale. 
So, we can conclude that (/) is constant, i.e., it does not depend on Q. This conclusion 
motivates the investigation for superior moments of /. 
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Similarly we can compute the second moment: 

/ /N tr ^ 2X 
(I") = ( IE* 



1=1 , , 

N t r,T c ,T e 



I N tr _ 9 N tr N tr 

(E* +EEVi 

\ 1=1 1=1 j=l 

N tr ,T c ,T e 



oo e -N N N tr / 

E -am - (M 7 / 

— N tr =0 ly tr- \ N ' Tc,T e 



+N tr (N tr - 1) {l)l T ) 
where (7J 7 ) = (7 7 ) = (7) 2 . 

\ ' J I T c ,T e \ l /T c ,T e \ J I T c ,T e \ I T c ,T e 

So, this amount can be calculated: 

7 2 \ - * f ( 1+e ' 9 ) 2 

c ' e — Q 



(2Q)(e«+l) 



(Q - e« + Qe« + 1) 



X2 Q-tanh(Q/2) 
2Q 



Because £~ r=0 ^^TTT = ^ and ^N tr =o N UN tr - 1)— = N 2 , thus we find that 

</ 2 > = ^|^[Q-tanh(g/2)] + (/) 2 

what gives var (/) = (/2) _ (/) 2 = ^_g!> [Q _ tanh(Q/2)] . And thus the relative error is 
described by 



-N e — N 



e(I) = ^> 



w 1 ■ 



Q-tanh(Q/2) 



(1) 



1/2 

O ( S2 ) \ (In 10 jVdee^ £)~ 1/2 r <3-tanh(Q/2) 1 1/2 
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= ^ ^ Q-tanh(Q/2) y 
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, , x 1/2 

\ S ) \ (In 10 N dec W L)~V 2 



where K = 2— s- i — 

\ {6) J (Pmax-Pmin) ' 

Thus by considering the variable q = (E? — ^i)jk^T being uniformly distributed, the 
relative error behaves as a simple logistic function of the band gap amplitude (Q). 

Aiming to propose a more realistic model then the supposition of the density of states 
following an uniform distribution, we propose an alternative model that considers the den- 
sity of states as being quadratically distributed. The theoretical quadratic function is a 
reasonable approximation to experimental results [8( and is simpler than the polynomial 
fit to those experimental data. For this, the theoretical u— shape modeled as a quadratic 
function defined in the interval [-Q, Q] is here proposed: 

Pr(g) = aq 2 + c. 

Supposing normalization J Q Q Pr(q)dq = 1, i.e., 2aQ 3 /3 + 2cQ = 1, we have 

, . , (1 - 2aQ 3 /3) 
Vr(q)=aq 2 + y - jQ^' (2) 

and from the existence condition c = (1 — 2aQ 3 /3) > and then a < 3/ {2Q 3 ) 

Following the same procedure used before (density of states uniformly-distributed) now 
applied to the case which the density of states follows a quadratic curve, we first calculate 
the first moment and then compute the average of current considering the random telegraph 
noise caused by all traps in the interface: 

(I) = (P/^S)P^ 

Q 



(l+e-i) an 2 I t 1 " 2 ^ 3 / 3 ) 
{S)N r 1 + ' 9 + 2 « 



(3) 

/ (l+coshg) ~dq 

Q 

This result must be obtained by using numerical methods for each a. For a — > 0, we have 
the expected (J) — > N £ ■ We tested many values a < 3/(2Q 3 ) and no meaningful deviates 
of (/) = N 2 were observed. Indeed just corrections of order O(10 -7 ) were verified, and 
there are no evidences to believe there might be large deviations for the current considering 
our theoretical u-shape curve. Our tests were performed using Q max = 10, and to perform 
it, just values a < 3/(2 -10 3 ) = 0.0015 were considered. We have employed the Simpsom rule 



(see for example 



10j) for the numeric integration of equation [31 The reader might think a 



7 



is very small and differences would be observed at higher values of a. In order to clarify this 
issue, we experimented a = 0.01 up to a = 0.1874 < 3/(2 ■ 2 3 ) = 0.1875 corresponding to 
Qmax = 2 and the same conclusion was obtained: (J) is constant under Q variations. Thus, 
the relationship between (J) and Q is independent of a. 

The most accurate model comes from the actual density of states, which can be measured 
from experimental data as in jg] and fitted by a polynomial. We have verified the differences 
in the MOSFET current fluctuations between modeling the band-g cip clS db uniform distribu- 
tion, a theoretical u-shape modeled by the simple quadratic function described by equation 
|2]and the actual experimental data fitted by a 8^-degree polynomial. 

The procedure of using the experimental density of states consists of performing a scan 
of the figures of the experimental density of states found in 8|, i.e., the density of traps 
having energy E T , corresponding to q — (E T — fi)/ksT, here denoted by f(q). 

In that work 3 u-shape curves are studied, fi(q) x q, /2(g) x q and, /3(g) x q, corresponding 
to 3 different oxides (in the nomenclature of ref. [8J], reoxidized: oxide I, nitrided: oxide II 
and, tceoxide: oxide III) defined in a symmetric band gap interval [—a, a]. The first step 
was to map this band gap to an interval [-Q, Q], what was done with the simple transfor- 
mation q' = Q 1 — ^ a ~ q ' . For our scanning n = 100 points, { (g;, /(<&)) }™ =1 , were col- 
lected from experimental figures and fitted by a 8-degree polynomial {gi{q') = Y^k=o@ k< i' k )- 
The degree of 8 was used because it presents a good visual fit for all u-shape curves. A 
normalization is required, and thus Prj(g') = gi(q')/ f Q Qgi(q')dq' = Ylk=i ^'kl' ' * ^ wriere 
P'k = N (2E£=o/Vu(2j + ^Y X Q 2J+1 )> f o r <?' e [~Q, Q]- And so a polynomial fit was 
performed for each value attributed to Q changing from Q m i n to Q ma x and after we have 
the experimental probability distribution describing the density of states. Next we per- 
form numerical integrations of the equations for the current and the relative error obtained 
in this manuscript, which are performed using Simpsom method with a very small step 
(h ~ 10~ 3 ). It supplies high precision for the integration and so for our estimates. These 
simple procedures are implemented in Fortran. 

Figure |5] shows the experimental u-shapes. The inside figure corresponds to an exam- 
ple of good polynomial fit performed for the one of the oxides (tceoxide). So we perform 
calculations to obtain (/). 

In figure |3] we can observe that (J) decays as function of Q which is not observed when 
the density of states is uniform or even when the curve is the theoretical u-shape. This 
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FIG. 2: Experimental density of states (U-shape curves) for 3 different oxides. The data was 
extracted from ref. 8]. The inside plot corresponds to a example of the good polynomial fit 
by a nine degree polynomial for the oxide III (tceoxide in the nomenclature of [§]). Just for a 
comparisom the continuous curve represents a theoretical curve for a = 0.09 for the same Q = 2. 



corroborates the idea that the asymmetry of the u-shape might be important to observe the 
current differences that appear when Q changes. However, it is important to notice that the 
average noise current decays at most 10% (for the tceoxide), for high values of Q. Since Q 
for Sj02 is approximately 2, in this case the decrease is nearly 5%. It is also important to 
analyze the relative error in the current by computing its higher moments. 
We then calculate the second moment for the theoretical u-shape density: 

Q 



Tc,T e 



(1 + e 



-q\2 



2 , (l-2aQ 3 /3) 



4(1 + coshg)^ 



dq 
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FIG. 3: Average current using the 3 experimental density of states shown in figure [2] as function of 
Q. The results show a decay which is not observed when the density of states is uniform or even 
when our theoretical u-shape density of states is considered. 



and the variance is 



var{I) = N(5 2 ) 



:i + e- 



-q\2 
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dq. 



We can calculate also in this case the relative error after a few algebraic manipulation we 
have 



e(I) = K- 
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dq 



I 



(1+e-i) 



2 (1-2^/3) 
a( f + 2Q 



(1+cosh q) 
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Integrating numerically the last equation we have 
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in the left plot of figure H] the behavior of e(I) as function of Q, for different value of 
a < 0.0015 showing that the relative error is sensitive for different and small (possible) 
values of a . For a suitable comparison, it is plotted in the same figure the case a = 0, 
which corresponds to the uniform distribution of density of states. The inside plot shows 
the same plot in a log-log scale. For higher values of a a power law behavior can be adopted 
to describe the e(I)x Q. The left plot of the same figure H] describes that the power law 
fit e(J) = aQ e considering Q > 2 is suitable simple fit for this interval. We find the best 
estimate = 0.143(1) and a = 1.013(1). 





FIG. 4: Left plot: Relative error as function of Q for a u-shape(quadratic) density of states. The 
inside plot shows a log-log plot for Q > 4 showing that higher values of a leads to a power law 
behavior of the relative error as function of Q. Right plot A power law fit of relative error as 
funcion of Q. It was used a = 0.00148 and Q > 2 for the plots. 

These results show that the relative error in the current is not universal since the relative 
error shows different behavior depending on the density of states. Finally we have performed 
the same calculations for the relative error considering the experimental u-shapes, as shown 
in figure [5j We observe a difference in the dependence of Q between the 3 different u- 
shapes. A similar behavior (logistic curve), predicted for the uniform case (see equation 
[1]), is observed using the experimental u-shapes, however there is a sensitive difference for 
higher values of Q between the two cases. Our theoretical approach for the u-shape supplies 
exactly an intermediate behavior between the uniform and experimental u-shape results, 
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showing it is suitable for modeling RTS noise in MOSFETs because it is simpler then the 
higher-order polynomial fit and more accurate than the uniform case! 




FIG. 5: Relative error as function of Q for the experimental u-shapes. 

Summarizing, we developed an original methodology in time domain to model the noise 
behavior due to traps at the Sj-Sj02 interface. Our results show that the relative error 
depends on the amplitude (Q) of the density of states at the interface. 

Our model for the average current and its relative error depends upon the model of 
the density of states adopted. When considering the density of states as having uniform 
or quadratic density of trap energies in the band gap we conclude that the relative error 
behaves as a logistic function in relation to Q. For the density of states modeled as an 
uniform distribution, a very simple formulation can still capture a logistic behavior similar 
to the experimental data. 

The experimental u-shapes show a decay in the average current not observed when the 
theoretical u-shape is considered. Such results represent an important advance in modelling 
RTS at device and circuit levels, and experiments can be performed to check these behaviors. 
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We believe that device designers can use the models proposed in this work for modeling the 
relative error as a function of the band-gap in upcoming technologies. 
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